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Abstract 

We have discussed the Tsallis entropy in finite A^-unit nonextensive systems, by using the mul- 
tivariate g-Gaussian probability distribution functions (PDFs) derived by the maximum entropy 
methods with the normal average and the g-average {q: the entropic index). The Tsallis entropy 
obtained by the g-average has an exponential dependence: Sq^"^ /N ~ e*-"*^"^^^ "^i ' for large 
N {i-q) > 0)- contrast, the Tsallis entropy obtained by the normal average is given by 
Sq'^^ /N ~ [l/{q — ^)N] for large N (q-i) ^ ^ dependences of the Tsallis entropy obtained 
by the q- and normal averages are generally quite different, although the both results are in fairly 
good agreement for \q — 1| ^ 1.0. The validity of the factorization approximation to PDFs which 
has been commonly adopted in the literature, has been examined. We have calculated correlations 
defined by Cm = {{Sxi 5xj)"^) — {{5xi)"^) {{6xj)'^) for i ^ j where 5xi = Xi — (xi), and the bracket 
(•) stands for the normal and (/-averages. The first-order correlation (m = 1) expresses the intrinsic 
correlation and higher-order correlations with m > 2 include nonextensivity-induced correlation, 
whose physical origin is elucidated in the superstatistics. 
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I. INTRODUCTION 



In the last decade, much attention has been paid to the nonextensive statistics since Tsal- 
hs proposed the so-called Tsallis entropy . The Tsallis entropy for A^-unit nonextensive 
systems is defined by 



^ 1 — q 



(1) 



where q is the entropic index, the Boltzmann constant {ks = 1 hereafter), pq^\x) 
denotes the A^-variate probability distribution function (PDF), x = {xi} (i = 1 to A^), and 
dx = YliLi dxi. The Tsallis entropy is a one-parameter generalization of the Boltzmann- 
Gibbs (BG) entropy, to which the Tsallis entropy reduces in the limit of g — ?■ 1.0, 



q{N) _ q{N) 



P. 



f)(^)lnpfH^) dx. 



(2) 



The BG entropy is additive because for two independent subsystems A and B where the 
PDF is assumed to be given by 



X2 



for xi E A, X2 G B, 



(3) 



we obtain 



si^\A + B) = si'\A) + s['\B). (4) 

In contrast, the Tsallis entropy is pseudoadditive (non- additive) [l-3, 6| because Eqs. ([I]) 
and (|3]) lead to 



Si'UA + B) = Sj^'\A) + Sj'\B) + il-q)S!'\A)S!^\B), 



(5) 



where the subscript of FA is attached for a later purpose. Similarly, when the PDF for 
A^-unit independent subsystems is assumed to be given in a factorized form, 

N 



(6) 



i=l 



we obtain the pseudoadditive Tsallis entropy, Sjf^Jj^, expressed by 



N 



In 



1 + (1 - <l)S, 



(N) 
q,FA 



5^ In [l + {l-q)S^'\^ 



(7) 



i=l 



We should note, however, that Eqs. ([5]) and ([7]) are not correct in the strict sense because 
the PDF derived by the maximum-entropy method (MEM) cannot be expressed by Eq. ([3]) 
or (E]), as will be shown shortly [Eqs. ( P7|) -( 15U]) ]. Indeed, for identical, independent systems, 
Eq. ([T]) with the use of exact multivariate PDFs yields 



C(2) 



'-'n — 'Jn FA ^ ^^rr 



q,FA 



(8) 
(9) 



where ASq' (= Sq"' — S'^"pj^) expresses a correction term. Equations ([HD and imply 



■,{N) l_ q{N) _ q{N) 

that Sq^^ does not satisfy the pseudoadditivity because the correction term of ASq'^^ is not 
negligibly small, as will be shown in this study. 

It is worthwhile to explain the nonextensivity which is a different concept from the non- 
additivity |7|. An entropy S^^^ for A^-unit system is said to be extensive if 



< lim 



< oo. 



(10) 



For independent (or short-ranged interaction) systems, the BG entropy is extensive: S^^^ 



BG 



NS^j^Q, while the Tsallis entropy is nonextensive. In contrast, the BG entropy is nonex- 
tensive for systems with long-ranged interactions. The Tsallis entropy may be extensive 



Gaussian PDFs with correlation 



io|. 



for correlated systems such as the probabilistic system [8|, some fermionic system and 



The PDF is commonly evaluated by the MEM for the Tsallis entropy with imposing some 
constraints. At the moment, there are four possible MEMs: (a) original method pj], (b) un- 
normalized method jj], (c) normalized method and (d) the optimal Lagrange multiplier 
(OLM) method [sl]. A comparison among the four MEMs is made in Ref. j^. In relation 
to the issue on the MEMs, two kinds of definitions have been considered for an expectation 
value of physical quantities: one is the normal average in the MEM (a) and the other is the q- 
average using the escort probability in the MEMs (c) and (d). Various arguments have been 



given that we should employ the g-average 6l.[ll-|l3|. Recently, however, it has been pointed 
out that for a small change of the PDF, thermodynamical averages obtained by the g-average 



are unstable whereas those obtained by the normal average are stable 14j-ll6|. In contrast, 
Ref. 171] has claimed that for the escort PDF, the Tsallis entropy and thermodynamical 
averages are robust. This issue on the stability (robustness) of thermodynamical averages 



as well as the Tsallis entropy is currently controversial 



18| 



In our previous papers [l9|,l20|, we discussed the Tsallis entropy and the generahzed Fisher 



information in nonextensive systems with and without spatial correlation, calculating the 
multivariate g-Gaussian PDFs. The purpose of the present paper is twofold:^(l) to make 
a comparison between the Tsallis entropies evaluated by the normal average [1] and the q- 
average 2|, |5|], and (2) to examine the validity of the FA and pseudoadditivity of the Tsallis 
entropy. One of the advantages of a use of the g-Gaussian PDF is that it is free from an 
open problem on ambiguity in defining the physical temperature in conformity with the 
zeroth law of thermodynamics 2l|- 25|. By using the derived PDFs, we may calculate the 
correlation induced by nonextensivity. Some related issues of the degree of freedom N on 
the PDF have been discussed in Refs. 



been calculated for classical ideal gas 



26 



29 



28|. The nonextensivity- induced correlation has 



31| and harmonic oscillator 30|. Our calculations 



will provide some insight to the current controversy on the normal versus g-averages and 
show the importance of effects which are not taken into account in the FA. 

The superstatistics is one of alternative approaches to the nonextensive statistics besides 



the MEM 



32 



- l34l | (for a recent review, see [35|). In the superstatistics, it is assumed that 



locally the equilibrium state of a given system is described by the Boltzmann-Gibbs statistics 



and its global properties may be expressed 
parameter {i.e., the inverse temperature) 



jy a superposition over the fluctuating intensive 



32]-[35|. The superstatistics has been adopted 



in many kinds of subjects such as hydrodynamic turbulence, cosmic ray and solar flares 
351] . The physical origin of the nonextensivity-induced correlation may be elucidated in the 
superstatistics. 

The paper is organized as follows. In Sec. II, we calculate the Tsallis entropy, by using 
multivariate PDFs for correlated nonextensive systems derived by the OLM-MEM with the 
g-average slQ, 2^. In Sec. Ill, we derive the multivariate PDF by the original MEM l| 
with the normal average in order to calculate the Tsallis entropy. In Sec. IV, we present some 
numerical calculations of the Tsallis entropy evaluated by the normal and g averages. In Sec. 
V we calculate correlation induced by nonextensivity, whose physical origin is elucidated in 
the superstatistics [32|, l33|. Sec. VI is devoted to our conclusion. 
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II. OLM-MEM WITH THE g- AVERAGE 



A. fl-Gaussian PDF 



We consider A^-unit nonextensive systems whose PDF, p^Hx), is derived with the use of 
the OLM-MEM js] for the TsaUis entropy giveribv Eq. ([I]) [l, 2|, 
given by (for details, see Appendix B of Ref. |20|) 



2|] . We impose four constraints 



1 



1 ^ 



i=l 
TV 



1=1 



N N 



s a 



MiV - 1) £ [{Xi - f^){xj - f^)]r 



(11) 
(12) 

(13) 
(14) 



Here /z, and s express the mean, variance, and degree of intrinsic correlation, respectively, 
and the bracket [-jg denotes the g- average over the escort PDF, P^q^\x), 



Pi''\x)Qix) dx, 



(15) 



with 



Pf)(a.) 



pf\x] 



I (pf )(..))^ dx, 



(16) 
(17) 



where Q{x) stands for an arbitrary function of x. 

The OLM-MEM with the constraints given by Eqs. ([II])-(II4]) leads to the PDF given by 



20 



p(f\x) 



AN) 



exp„ 



<^(x) 



(18) 
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where 



<l>(a;) 



ao 
ai 



N N 

y^[Qo + ai(l - - - yu), 

i=l i=l 

[1 + (iV - 2)^] 
il-s)[l + {N-l)sy 
s 



(l-.)[l + (iV-l).]^ 



(19) 

(20) 
(21) 



(N) \2^r4^ 
[ 9-1 



N/2 iv^ 



1-g 



r(T^+f+i) 



for g > 1, 
for q = I, 
for g < 1, 



{(l-s)^-i[l + (iV-l)s]}V2, 



« 2 



(22) 

(23) 
(24) 



Here B{x,y) and r(2) denote the beta and gamma functions, respectively, and expg{x) 
expresses the g-exponential function defined by 

exp^(x) = [I + (I - q)x]f'-'^\ (25) 

with [a;]_|_ = max(x, 0), which becomes expg(x) — )■ for q — )■ 1.0. The entropic index q may 
take a value, 



0<q<l + — = qu, 



(26) 



because p[^^ (x) given by Eq. (fT8l) has the probability properties with z/g^-* > for g < g^/ 
and because the Tsallis entropy is stable for g > Evaluations of Zq^^ and [Q]q with 
the use of the exact approach |40|, |4l| are discussed in Appendix A.l. 



In the absence of the intrinsic correlation (s = 0), Eq. fITS]) reduces to 19| 



7{N) 



exp„ 



N 



2z.rv^ 



i=l 



(27) 



On the other hand, the PDF in the FA is given by 



N 



Up 



i=l 



1 ^ 



exp„ 



i=l 



24'^^ 



(28) 
(29) 



In Eqs. fl27|) and (|29|1 . Zq is given by Eq. f l22|) with rs = 1.0. From a comparison 
between Eqs. (l27|l and (129|1 . it is evident that 

^ pr^^^(^)' (30) 

except for g = 1.0 or = 1. 



B. Tsallis entropy 

Substituting the PDF given by Eqs. ( !T8|) - (l2^ to Eq. ([1]), we obtain the TsaUis entropy 



given by 



19 



20| 



1 - 



g-1 

with 



= (31) 



cr = -r[^r]'"- (32) 



The s-dependence of the Tsalhs entropy was previously discussed (see Fig. 1 of Ref. |20|). 
With increasing s, the Tsalhs entropy is decreased as given by 

SW{s) = 5f )(0) - ~ ^^^^"^^ for \s\ « 2/v/iV(iV-l) . (33) 

We pay hereafter our attention to identical, independent systems with s = 0. The Tsallis 
entropy for q = 1.0 (BG entropy) is given by 

S-;"^) = NS['\ (34) 

with 

) = [ln(27ra') + 1] = 1.418 for = 1.0. (35) 



Employing Eq.( !29|) . we obtain the Tsallis entropy in the FA given by 

9-1 



Sl% = i^^, (36) 

N 



iVSW + ^(^-i)(i--') (5ay + ... (38) 
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In particular for N = 2, Eq. (138|) becomes 



1-fc 



,(1)a2 



^•^^ q - 1 

By using the formula: 



S 'L = = 2^i^) + (1 - <lM^r. (39) 



, , T{z + a)\ a(l - a) . , , 
In — ^— — - ~ ^ for hz oo, (40) 



we obtain the entropy for |g — 1| ^ 1.0 given by 



(41) 



sr - sfA - A's!" - (9 - 1) - ^ 

The (g — 1) term in Eq. fHT]) includes 0{N'^) contributions showing the nonadditivity. 
For large iV, we obtain the Tsallis entropy expressed by 

^ (^^j (42) 

employing the relation: lnr(2;) ~ (z — 1/2) In 2 — 2; + (1/2) ln27r + ■■ for | z\ — }■ 00. Sq ^ given 
by Eq. ( H2l) leads to a very large figure. In the case of g = 0.5 and = 1.0, for example, 
4^V^ is 5.99 X 10^, 3.09 x lO^o and 4.11 x 10^°^ for = 10, 100 and 1000, respectively. 
It yields an astronomical figure for Avogadro's number of A^ = 6.022 x 10^^. 

III. ORIGINAL MEM WITH THE NORMAL AVERAGE 
A. g-Gaussian PDF 

In preceding Sec. II, we have made calculations by using the OLM-MEM with the c- 
average jsl. In this section, we employ the original MEM with the normal average 
imposing the constraints, 

''p[^^\x)dx, (44) 
1 ^ 

(45) 



/i 



N 

i=l 



1 ^ 

-5^((x.-/i)\, (46) 



N 

i=l 



^ N N 

NiN -I) ^ 5Z (47) 



where the bracket (■)g denotes the normal average (relevant quantities being expressed with 
tilde hereafter), 



The original MEM yields 



p^f\x)Q{x) dx. 



with 



P[^\x) 



Exp, 



^ix) 



(<?-!) 

rf)(27ra2)^/2 



^9-1 2 



1 N/2 pf^.iv^ 



. 1-9 2 



^ for g > 1, 
for g = 1, 
for q < I, 



(48) 



(49) 



(50) 



2 I g 



1 = U- 



(N) 
1/9- 



Here $(a3), rs^-* and z/g^^"" are given by Ens. f[T9l) . fl23|) and flM|) . respectively, and Exp (x) is 



(51) 



the g-exponential function defined by [l 

Expg(x) = 







1 + 









1/(9-1) 



(52) 



which is different from expq(x) given by Eq. ( I25p . The two g-exponential functions, exp^(a;) 
and Expg{x), have the relation jl5 |: 

expg(x) = Exp2_g((2-g)x), Expg(x) = exp2_g(a;/g). (53) 



The condition of ujj^'' > implies that a conceivable g value is 

2 1 

N + 2 qu 



(54) 



where qu is given by Eq. (!26|) . The PDF given by Eqs. (H9l) - (!5T!) is flat-tailed for g^, < g < 1 
and compact support for g > 1, which is in contrast with the PDF given by Eq. fl271) obtained 
by the OLM-MEM. Evaluations of Zg^^ and {Q)q with the use of the exact approach 40, 4l| 
are discussed in Appendix A. 2. 
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B. Tsallis entropy 



By using the PDF given by Eqs. (H9|) - (l5T|) . we obtain the Tsalhs entropy given by 



n{N) 



1-5, 



q-l 



with 



1-q 



(55) 

(56) 
(57) 



The Tsalhs entropy in the FA is given by 



'~'q,FA 



i-[5; 



;;(l)liV 



q-l 



For |g — 1| ^ 1.0, the entropy is given by 



~ S. 



(N) 
q,FA 



NS['^-{q-l] 



+ 



which is the same as Eq. (HTl) except for the 0{N) term in the bracket. 
For large N, the Tsallis entropy is expressed by 



r.(7V) 



(g-l) 



2q e-(^"i)^^i 
N{q-1) 



(1) 



s. 



(N) 
q,FA 



^ ^ for iV » ^ > 0, 



q-l 



q-l 



which lead to S[^^ 



g-dependences of Sq^"* and S^^^'p^^ in the normal average are quite different from those of 5"^ 
[Eq. (|42])] and ^^^^^ [Eq. ^] in the g-average. 



and A5f ) 



(58) 



(59) 



(60) 
(61) 



for large (g — 1)A^. We note that the A^- and 



(JV) 



IV. NUMERICAL CALCULATIONS 

We show some model calculations in which we hereafter set /x = 0.0 and o"^ = 1.0. 
Equations (l26l) and (|5^ show that the results obtained by the q- and normal averages are 
valid forO<g< 1 + 2/N = qu and q > 1 — 2/(A^ + 2) = g^, respectively (see the inset 
of Fig. [1]). Although for large A^|g — 1|, the A^ dependence of the Tsallis entropy of the 
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g-average [Eq. ( H2|) ] is quite different from that of Sg^^ of tlie normal average [Eq. ( 160|) ]. 
botfi entropies are in fairly good agreement for < q < qu. Figure 1 shows the Tsalhs 
entropy for = 10 as a function of q calculated by the g-average (solid curve), normal 
average (chian curve) and by the FA in the normal average (dashed curve). The differences 
between them are small except for q ^ qi where Sq^^ divergently increases. 

Figures 2(a) and 2(b) show three-dimensional plots of Sq^^ /N and ASq^^ / Sg^\ respec- 
tively, calculated by the g-average as functions of q and A^, where ASq^^ = Sq^^ — S^^Jj^. 
We note in Fig. 2(a) that Sq^^ is exponentially increased with increasing (1 — q) and/or 
A^. Figure 2(b) shows that ASg^^ / Sg^^ — )■ 1.0 for large (1 — q)N, where the Tsallis entropy 
IS given by ^ ASi^K This implies that the FA is not a good approximation and the 
pseudoadditivity is violated. 

Figures 3(a) and 3(b) show three-dimensional plots of Sg^^/N and ASg^^ / Sg^\ respec- 
tively, calculated by the normal average as functions of q and A^, where ASq^^ = Sq^^ — S^^J^. 
It is shown that Sq^^ is decreased with increasing (g — 1) and/or A^, where ASq^^ / Sq^^ is 
considerably decreased. The maximum value of ASq^^/Sq^'^ is about 0.002, which means 
that the pseudoadditivity of Sg^^ is approximately satisfied with an accuracy better than 
0.998, and then the FA is a good approximation. 



V. DISCUSSION 
A. Correlations 



1. q- average 

We will calculate the mth-order correlation for i ^ j defined by 

Cm = [{Sx, 5x,r]q - [i5x,r]g [{5xir]q, (62) 

which is evaluated by the g-average {6xi = Xi — n). With the use of the PDF given by 
Eqs. (IT8l)-( IMl) . we obtain the first- and second-order correlations given by (for details, see 
Appendix A.l) 

Ci = [5xi Sxj]q = cr^s, (63) 

C2 = C2s + C2„, (64) 
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with 



2[(A^ + 2) - A^g] (7^2 
[{N + 4)-{N + 2)q] ' 
2{q - 1) a' 



C2S 

C2n 



for < g < 1 



Af+2' 



(65) 
(66) 



[(AT + 4) - {N + 2)q] 

where we adopt [6xi]q = 0. We note that Ci and C2S arise from intrinsic correlation s, and 
that expresses correlation induced by nonextensivity which vanishes for g = 1.0 and 
which approaches —2a'^/N as N ^ 00. In particular for N = 2, C2S and C2n are given by 

2(2 - q)a\s^ 



C2S 

C2n 



(3 - 2g) 
(g - 1) 



for < g < |, 



(3 - 2g) 

where > for < g < 1.5 while C2n < for < g < 1 and C2n > for 1.0 < g < 1.5 
On the contrary, the PDF in the FA given by Eq. (129|) yields 



(67) 
(68) 



FA 



'^2s — 



FA 
2n 



0. 



(69) 



By using Eq. (127|) . we obtain the mth-order correlation for arbitrary m with s = 0, 



Cr, 



Al 



Al 



(9-1) 



r(^-f->») 
r(^-f) 



r(T^+f+i) 
r(T^+f+™+i) 



\q-l 2 2 J 

r(^-f) 



r(i^+f+i) 
r(i^+f+f+i: 



for 1.0 < g < g™, 
for g < 1.0, 



(70) 



where 



Ar, 



-jy-^ for even m. 







1 + 



for odd m. 



(71) 



(72) 



A^ + 2(m- 1) 

The nonextensivity yields higher-order correlation of Cm with m > 2 for g 7^ 1.0 in indepen- 
dent nonextensive systems where s = Ci = 0.0. 

Figure 4 shows C2 as functions of g and s for = 2. At the origin of (g, s) = (0.0, 0.0), 
we obtain C2 = —0.33. The parabohc-like line starting from (g, s) = (1.0,0.0) expresses the 
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intersection of C2 with the zero-level plane. With increasing q and/or s, C2 is monotonously 
increased. As q approaches 1.5, C2 is divergently increased. 

Figure 5(a), 5(b) and 5(c) show Cm with m = 2, 4 and 6, respectively, as functions of q 
and A^. Magnitudes of Cm are significantly increased with increasing m. Cq shows peculiar 
q and dependences. 

2. Normal average 



We adopt the normal average to evaluate the mth-order correlation for j 7^ j defined by 



(73) 



With the use of the PDF given by Eqs. (149|) -( 15T|) . the first- and second-order correlations 
are given by (for details, see Appendix A. 2) 



Ci = {5Xi SXj)q = cr^s, 
6*2 = C2S + C2n, 



with 



2[(A^ + 2)g- A^] a^s^ 
{N + 4)q-{N + 2) 
2(1 -g) a' 



C2S 

C2n 



for g > 1 



(AT + 4)g - (AT + 2) ' ^ ^+4' 

where {Sxi)q = 0. In the case of A^ = 2, Eqs. (!76|) and (1771) become 

2(2g - 1) 



C2S 

C2n 



3q-2 

^ — — for g > |. 

3g - 2 ^3 



Cn 



At 



(g-i) 



, q-1 ^ 2 



+1 



Al 



(1-9) 



r(i^+f+-+i) 



\l-q 2 

r(T^-f) 



\ 9-1 ' 2 ' / 
r(i^+f+T+i) 



(74) 
(75) 



(76) 
(77) 



(78) 
(79) 



By using Eq. fl49|) . we obtain correlation of Cm for arbitrary m with s = 0, 

for g > 1.0, 



'^L JV _ m'l 

.1-9 2 2 > 



, 1-9 2 



for g < 1.0, 



(80) 
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where Am is given by Eq. (17T]) . When comparing Eqs. (^^, (176ll . (177|) and (IHOj) with 
Eqs. fl63|) . fl65|) . fl66|) and flTOj) . respectively, we note that Cm and Cm have the reciprocal 
symmetry: q 1/q. Then q dependences of Cm are given by those of Cm in Figs. 4 and 5 
if we read g — )■ 1/g. 



B. PDF in the superstatistics 



The physical origin of the nonextensivity-induced correlation is easily understood in the 



32- 



superstatistics 
composite systems 
given by [2(1] 



34 



It is straightforward to apply a concept of the superstatistics to 



10|. We consider the A^-unit Langevin model subjected to additive noise 



dxi 



-Xxi + 72/} ii{t) + 1 for i = 1 to A^, 



(81) 



where A denotes the relaxation rate, ^j(t) the white Gaussian noise with the intensity D, 
and / an external input. The PDF of tt^'^\x) for the system is given by 



N 



TT 



(N) 



(82) 



i=l 



where the univariate PDF of iT^^\xi) obeys the Fokker-Planck equation. 



-^[{Xxi - I)n^'\xi, t)] + D^TT^'^ixi, t). 



92 



dt Sx, 
The stationary PDF of 7c^^\xi) is given by 

1 



dx. 



(83) 



V2 



: exp 



TTCr^ 



2^2 



(84) 



with 



fi = I/X, a^ = D/X. 



(85) 



After the concept in the superstatistics 32h35| . we assume that a model parameter of (3 
= rA/W. fluctuates, and that its distribution is expressed by the x^.^jig^rji^ution with rank 



n 



32 



33| 



r(n/2) V2/3, 



n 



n/2 



__ j ^n/2-lg-n/3/2/3o^ 



(86) 
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where T{x) is the gamma function. Average and variance of /3 are given by = /3o and 
((/^^)/3 ~ f^o)//^o — '^I'^i respectively. Taking the average of -k'^^^x) over we obtain the 



'/3 

stationary PDF given by 



19| 



oc 



1 + ^5^(0:.-/.)^ 



N 



1=1 



'{N+n)/2 



which is rewritten as 



with 



7(N) 



exp„ 



/3o 



N 



2u, 



(TV) 



i=l 



{N) 1 TV 



(g-l)/9o 
2 

^ (A^ + n) 



i=l 



1 1 



2' g- 1 



(87) 
(88) 

(89) 

(90) 
(91) 



Uq^^ being given by Eq. (j24l) . In the hmit of n — )• oo (g — )■ 1.0) where /(/3) — )■ — /^o), the 



PDF reduces to the multivariate Gaussian distribution given by 



pf\x) 



1 ^ 



(92) 



which agrees with Eqs. ([82]) and ([84]) for /3o = A/L* = l/o-^. 

We note that the PDF given by Eq. (189)) is equivalent to that given by Eq. (127)) derived 
by the MEM when we read Pq = 1/cr^, although the former is valid for l<g<l + 2/ (A^ + 1) 
while the latter for < q < 1 + 2/N. The nonextensivity-induced correlation arises from the 
common fluctuating field of /3, because pq^\x) ^ WiVq^i^^i) for g 7^ 1.0 despite 7r*^^^(cc) = 

Alternatively we may rewrite the PDF given by Eq. ( 188)) as to conform to the PDF in 
the normal average for s = [Eq. (1^9)) ]. 



pf^[x) oc Exp = 



TV 







r,(^-2) 



with 



2i>. 



1 - 



4=1 



(93) 



(94) 



AT + n' 

where z/g^'' is given by Eq. (15T)) . The PDF given by Eqs. (193)) and (IM)) is valid for 
l-2/(A^ + n) < g< 1. 
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TABLE I: A comparison between the Tsallis entropies obtained by the q- and normal averages 



[Sj^^ = (l/2){ln(27r(j2) + 1}, = 1 + 2/N and qi = I - 2/{N + 2) = 1/qu]. 



method 


range 


Sg (for |g - 1| < 1.0) 


Sq (for iV » 1) 


g- aver age 


0<q<qu 






2 


normal average 


q> QL 






(g-1) Af(g-l) "2 J 



VI. CONCLUDING REMARKS 

Table 1 shows a comparison between the q- and normal averages, which have been ob- 
tained with the use of exact A^-variate g-Gaussian PDFs. Our calculation has shown the 
foUowings: 

(i) the MEMs with the q- and normal averages are valid for 0<g<g;7 = l + 2/A^ and 
q > qL = 1 ~ 2 / {N + 2) , respectively (see the inset of Fig. 1), 

(ii) the Tsallis entropy of Sg^'^ obtained by the g-average [sj has an exponential A^ dependence 
for large A^ where the FA is not a good approximation and the pseudoadditivity does not 
hold because ASjj^^ (= Sq'^'^ — S^^J^) is considerably large, 

(iii) the A^ dependence of Sq'^^ for large A^ obtained by the normal average is quite different 
from that of Sq^^ obtained by the g-average, and the FA in the normal average becomes a 
good approximation because ASg^^ (= Sq^^ — S^^Jj^) is very small, 

(iv) the q- and normal averages yield nearly the same results for qi < 1 < Qu where the two 
averages are valid, 

(v) the nonextensivity-induced correlation is realized in higher-order correlations {Cm and 
Cm) with m > 2 while the first-order correlation expresses the intrinsic correlation, and 

(vi) the nonextensivity-induced correlation is elucidated as arising from common fluctuating 



field introduced in the superstatistics j32l-l34|. 



Items (i)-(iv) clarify the difference and similarity between the q- and normal averages. Items 
(ii) and (iii) may suggest that the normal average is more favorable than the g-average, which 
is consistent with Refs. We should be careful in employing the FA with the 

g-average which is shown not to be a good approximating method except for |g — 1| ^ 1.0 
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although it has been widely adopted in nonextensive classical statistics 37|]. The situation is 



;he same also in the FA for nonextensive quantum systems as recently pointed out in Refs. 



38 
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Appendix: A. Evaluations of averages by the exact approach 



1. The g-average 



We first discuss an evaluation of the g-average given by 
ZrW./|l-(l-,)a.*(.)]ri..., 

Qj"l(a) = lQ{x)],= ( ' ) / Q(x) [1 - (1 - <,)a Hx)]^ dx, 

Vq Zq J 

by using the exact expressions for the gamma function 

1 



38 



40 



4l|: 



y 



for s > 0, 



-^r(s + 1) / {-r"^)e-y' dt for s > 0. 

ZTT 



(Al) 
(A2) 



(A3) 
(A4) 



Here a = l/{2vf\'^), $ {x) is given by Eq. f|T9l) . Q{x) denotes an arbitrary function of x, 
C the Hankel path in the complex plane, and Eq. (J32l) is employed. We obtain [381. kol. |41| 



Zf)(a) 



-tVy /o^ M^"^e-"Zj^^ [(g - l)au] du for q > 1.0, 

[q-l) 

+ l) /^(-t)-T^-ie-*zf )[-(l - q)at] dt for q < 1.0, 



(A5) 
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WT^WYh'^ JT '^"'^ ""^i^^iQ - 1)H Qi^\i(l - 1)H du for q > 1.0, 



^1 ^1 '-yq^i 



X qS^^[-(1 dt 



for g < 1.0, 
(A6) 



where 



Z\ '[a) J 



(A7) 
(A8) 



Thus we may evaluate the q average of Q{x) from its average over the Gaussian PDF. 
For example, simple calculations lead to 



with 



{{xi - = — 



1 

2a' 

s 



{{Xi - fl){Xj - = — 



2a 

A 



for i j 
for i ^ j. 



Q,m,/2 ' 

Al 



((x. - /i)"(a;, - iiT)i = ^ for ^ ^ J, s = 



a' 



(A9) 
(AlO) 

(All) 

(A12) 

(A13) 



for even m, 



for odd m. 



(A14) 



Employing Eqs. (jMl), dM]) and (lA9l)-( lAT3|l . we obtain 



(A15) 



for i 7^ j. 



(A16) 



2, (iV + 2-iVg)(l + 2s2)a^ 



(AT + 4) - (Ar + 2)g 



for^^j, (A17) 
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< 










1-g 



r(T^+f+f+i) 



for g > 1, 
for g < 1, 



(A18) 



and for s = and i ^ j, 
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2.rv2^^(^-f-„^) 



g-1 



1-g 



, l-q ' 2 



r(T^+f+-+i) 



for g > 1, 
for g < 1, 



(A19) 



which yield Eqs. (1^. (IBS]). (IBBj) and (170]). 



2. The normal average 



Next we discuss an evaluation of the normal average given by 
Z^g^\a) = J [l-{q-l)a^ix)]^ dx, 

Q{x) [1 - (g - 1) a dx, 



where a = 1 / {2qv^^ a"^) . By using Eqs. (1A3I1 and ( 1A4I) . we obtain 



(TV) 



(A20) 
(A21) 



+ 1) /^(-t)-^-ie-* )[-(g - l)at] dt for g > 1.0, 
r[i/(\-g)] M^~^e-" Z[^\{1 - q)au] du for g < 1.0, 



(A22) 



'^^'^^Pl^'^ ^[-(g - l)at] gi[-(g - l)at] t/t for g > 1.0, 



1 



^ zrnm-,-] io'^'^^'e'" - Qiia - qM du for g < 1.0, 



(A23) 



where z[^''{a) and (^^^^''(a) are given by Eqs. flA7l) and flA8l) . respectively. 



l(JV) 
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Employing Eqs. f R9|) - fRT3|) . flX22|) and fE23l) . we obtain 



a 



[{N + 2)q - N]{1 + 2s^) 
{N + A)q-{N + 2) 



for i ^ j, 
for i ^ J, 



(A24) 
(A25) 
(A26) 



q-1 
l-<? 



W-i^ 2 for > 1 

m/2 p/'^L iv_m^ 

'-\l~q 2 2y 
r(l^-f) 



for q < 1, 



(A27) 



and for s = and i ^ j, 



/12 

m 



(7-1 



, 9-1 ' 2 



+ 1 



1-9 



for g > 1, 



for g < 1, 



, 1-9 2 



(A28) 



which lead to Eqs. (|7i|), (j76|), (1771) and (|80|). 
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FIG. 1: (Color online) The g-dependence of Sq /N calculated by the normal average (N-av: solid 
curve) and g-average (q-av.: chain curve) and in the FA by the normal average (dashed curve) 
with = 10 for which qi = 0.833 and qjj = 1.20. The inset shows the validity range of the two 
averages in the q-N space: the g-average and normal average are valid for < q < qjj and qi < q, 
respectively, and for qL < Q < Qu the both averages are valid: the dashed line expresses = 10. 

FIG. 2: (Color online) (a) Sq^'^ /N and (b) ASq'^'^ / Sq^'^ of the Tsallis entropy calculated by the q 
average as functions of q and A^ (s = 0): note that the ordinate of (a) is in the logarithmic scale. 
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